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The chiral magnetic effect is the generation of electric current of quarks along external magnetic 
field in the background of topologically nontrivial gluon fields. There is a recent evidence that this 
effect is observed by the STAR Collaboration in heavy ion collisions at RHIC. In our paper we 
study qualitative signatures of the chiral magnetic effect using quenched lattice simulations. We 
find indications that the electric current is indeed enhanced in the direction of the magnetic field 
both in equilibrium configurations of the quantum gluon fields and in a smooth gluon background 
with nonzero topological charge. In the confinement phase the magnetic field enhances the local 
fluctuations of both the electric charge and chiral charge densities. In the deconfinement phase the 
effects of the magnetic field become smaller, possibly due to thermal screening. Using a simple 
model of a fireball we obtain a good agreement between our data and experimental results of STAR 
Collaboration. 
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I. INTRODUCTION 

Properties of strongly interacting matter in hadron- 
scale magnetic fields have attracted a lot of attention 
recently. The interest is motivated by feasibility to cre- 
ate such strong fields in ongoing heavy-ion experiments 
at Relativistic Heavy Ion Collider (RHIC), and at fu- 
ture experiments at Facility for Antiproton and Ion Re- 
search at GSI, and ALICE experiment at LHC. Non- 
central heavy ion collisions may create magnetic fields 
due to relative motion of electrically charged ions and 
the products of the collision. For typical RHIC parame- 
ters the strength of the magnetic field at the center of a 
Au-Au collision may be estimated to be of the order of 
eB ~ (10 - 100 MeV) 2 [1] at first moments (r - 1 fm/c) 
of the collision. The magnetic field created at RHIC 
is not only many orders of magnitude stronger than 
any other experiment may reach, it also leaves behind 
even magnetized neutron stars, magnetars, with their 
magnetic fields of the order of 10 10 Tesla (equivalent to 
eB ps 4 MeV 2 ). Strong magnetic field can significantly 
modify the properties of nuclear matter affecting the 
spectrum of the hadronic states Q, shifting the posi- 
tion and changing the order of the phase transition from 
hadronic matter to quark-gluon plasma [3[ . 

The magnetic fields lead to other unusual effects due 
to nontrivial topological structure of QCD [l], 0, [H, 
In this paper we concentrate on a particular realiza- 
tion of the Chiral Magnetic Effect (CME), which gener- 
ates a spatial separation of positive and negative electric 
charges along the direction of the magnetic field on non- 
trivial topological backgrounds of gluons [2, H[. In the 
noncentral collisions the magnetic field is perpendicular 
to the reaction plane, therefore a charge separation can 
be observed in heavy-ion collisions as a non-statistical 
asymmetry in the number of positively and negatively 



charged particles emitted on different sides of the reac- 
tion plane @, 0] • This effect was called "event-by-event 
V- and CP-violation" , because this asymmetry implies 
that the difference of the numbers of quarks of differ- 
ent chiralities was created in a particular event due to 
quantum fluctuations of the topological charge. However, 
since there is no CV - violation in strong interactions in 
the usual sense, this effect vanishes after averaging over 
all events. The only way to observe it experimentally 
is by studying CP-even correlations of reaction products 
Q. The magnetic field grows with the impact parame- 
ter of the collision, therefore this asymmetry should de- 
pend strongly on the centrality of the collision. There 
are preliminary indications that this effect has been in- 
deed observed by the STAR collaboration at RHIC 0,(1]. 
It should be stressed that the observed effect is the in- 
dication that there are fluctuations of chirality in QCD 
vacuum. A similar well-known effect is the large mass of 
the rf meson fiof . However, the CME allows one to 
extract more information on the fluctuations of chirality 
and topological charge, since in this case there is an ad- 
ditional parameter - the strength of the magnetic field - 
which can be varied in experiments. 

The physical mechanism behind the CME is as fol- 
lows [H, 0j- At the energy scales of the collisions, the 
light u and d quarks may be considered as approximately 
massless particles. The massless fermions are character- 
ized by right- or left-handed helicity (i.e., positive or neg- 
ative projection of fermion's spin on the momentum, re- 
spectively), so that the magnetic moment of (anti)quark 
is always collinear to its momentum. The strong enough 
external magnetic field makes the magnetic moment to 
be parallel to the direction of the field so that the mo- 
tion of the quarks is essentially collinear to the magnetic 
field. The helicity is related to chirality: the quarks have 
the same helicity and chirality while for the antiquarks 
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the helicity is opposite to chirality. In the equilibrium 
ensembles the total chirality is vanishing [TTf| , so that the 
total electric current along the magnetic field is zero (the 
left-handed and right-handed quarks move in opposite di- 
rections and compensate each other). However, if there 
is a difference between the (anti)quarks with left- and 
right-handed chirality then the balance is broken and the 
total electric current a long the magnetic field becomes 
nonzero 0, 0, i, [H [H [3] . 

Thus, the CME arises due to the fluctuations in the 
numbers of left- and right-handed quarks in a vacuum 
or in thermal state of non-Abelian gauge theories. In 
QCD the fluctuation of the chirality number are provided 
by the topology-changing transitions between different 
vacua. Indeed, the gluon fields are characterized by the 
topological winding number [IH [l6| 

Q = ^Jd 4 xG^(x)G a ^(x), QeZ, (1) 

where G^ is the gluon field strength tensor and the dual 
tensor is G» v = (1/2) e^ va ^G a ^. The configurations 
with Q ^ change the chirality of the ensemble, i.e. 
the difference between the number of left-handed quarks 
(Nl) and right-handed (Nr) quarks by an integer, 

A(N L - N R ) = NfQ, (2) 

where Nf is the number of flavors. Therefore the electric 
current along the strong enough magnetic field should 
be sensitive to the topological fluctuations of the gluon 
ensembles, providing us with possibility to study in the 
direct way the topological structure of the QCD thermal 
states. 

The CME is expected to work effectively in the chirally 
restored phase [3,13]: the chiral condensate, which breaks 
the chiral symmetry, provides a direct coupling between 
left-handed quarks to the right-handed ones thus leading 
to an efficient mechanism to suppress any asymmetry be- 
tween them. The suppression happens because the sys- 
tem prefers to decrease total chirality as the states with 
nonzero chirality have higher energy due to the Fermi- 
principle [ll| . Therefore the CME should work effectively 
in the high temperature case, where the chiral condensate 
is suppressed and the chirality fluctuations are strong. 

While there are many phenomenological estimates of 
the strength of the CME in QCD, a detailed analytical 
calculation is hardly possible due to strong nonperturba- 
tivc effects. In this paper we report on our results of the 
numerical investigation of the CME in quenched SU (2) 
lattice gauge theory both in high and low temperature 
phases of the theory (Section [TT|. The chiral properties 
of the QCD, such as chiral condensate [ijj and the chiral 
magnetization [181 ] , in the external magnetic field can rea- 
sonably be described by the quenched lattice theory. Our 
study involves various approximations: we use quenched 
lattice gauge theory with two colors, the fermionic propa- 
gators are truncated, and the theory is studied in a finite 
volume. Thus our results should only serve as a quali- 
tative indication of the existence of the CME, while the 



quantitative features may be accessible only in simula- 
tions with dynamical quarks. Using our numerical results 
we nevertheless make a prediction for the charge asym- 
metry and fit it to the experimental data in Section IIIII 
Our conclusions are summarized in the last Section. 



II. LATTICE SIMULATIONS 
A. Numerical setup 

We work in SU (2) lattice gauge theory with tadpole- 
improved Wilson-Symanzik action (see, e.g., expres- 
sion (1) in Ref. [191]). In our zero-temperature simula- 
tions we have used 14 4 and 16 4 lattices at lattice spacing 
a = 0.103 fm and a = 0.103 fm, a — 0.089 fm, respec- 
tively. The latter lattices were used to control the effects 
of finite volume and finite lattice spacings, which turn 
out to be much less than the statistical errors. For simu- 
lations at finite temperature we have used 16 3 x 6 lattices 
with spacings a = 0.128 fm and a = 0.095 fm, which cor- 
respond to the temperatures T — 256 MeV = 0.82T C and 
T = 350 MeV = 1.12T C , respectively. The critical tem- 
perature in SU(2) gauge theory is T c = 313.(3) MeV [2(J. 
As we explain below, the minimal nonzero value of mag- 
netic field is rather large in our simulations, <7-B m i n ~ 
(350 MeV) 2 ). 

In all current lattice simulations with background mag- 
netic field only valence quarks interact with electromag- 
netic field 0, [TtI HH [Hj], because the inclusion of dy- 
namical sea quarks into simulations make the problem 
significantly difficult. On the other hand the characteris- 
tic properties of the CME may be studied in the quenched 
approximation, in which the effects of the virtual quarks 
on the gluon fields are neglected. 

The quenched world has the same qualitative features 
as the unquenched system: both approaches exhibit color 
confinement and chiral symmetry breaking. There is, 
however, an important difference: due to the Crewther 
theorem [23| fluctuations of the net topological charge in 
the theory with light virtual quarks should be drastically 
suppressed compared to the quenched vacuum. However, 
as we have checked numerically, at least in the quenched 
approximation the CME is mostly due to local fluctu- 
ations of the density of topological charge, rather than 
due to the global topological charge of gauge field con- 
figurations. The number of colors is also inessential for 
the existence of the CME. Thus we proceed with the 
quenched simulations of the SU (2) gauge theory on the 
lattice. 

For valence quarks we use Neuberger's overlap Dirac 
operator [24[ . The expectation values of fermionic fields 
are first computed in terms of massive Dirac propagator 
in a fixed configuration of gauge fields, according to the 
Wick theorem. For example, for the expectation value 
involving four fermionic fields in background of a fixed 
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jauge field A, one has: 
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where V is the massless Dirac operator on the lattice, 



(4) 



and Oi, i = 1,2 are some operators in spinor and color 
space. These operators are treated as matrices in spinor 
and color spaces. Consequently, the traces in Eq. ((3J) are 
taken over both spinor and color indices. 

The overlap operator is a vital tool for investigation 
of the CME because this effect is realized for the chiral 
quarks, while the overlap operator is the lattice Dirac op- 
erator which enjoys the chiral symmetry by the construc- 
tion. On the contrary, if the quarks are massive then the 
chiral magnetic effect is overshadowed by the flips of the 
quarks' spirality. Moreover, in the CME the near-zero 
Euclidean eigenmodes play an essential role. The over- 
lap operator gently treats low-lying eigenmodes making 
our results more reliable. On the other hand, this oper- 
ator is more complicated compared to other versions of 
the lattice Dirac operators. As a result, the simulations 
are more power-consuming and our data has larger error 
bars. 

The Dirac propagator 1/(1? + m) can be represented 
in terms of the eigenvalues Xk and the eigenmodes ipk of 
the Dirac operator, 



Vip k = iX k tp k 



as follows: 



7 — 0,2/) = V 

4_ rm c * 



V + m 



fk (x) ijjj (y) 
i\k + m 



(5) 



(6) 



Then Eq. (|3|) can be reformulated as follows: 



= (k\Oi\k)(p\0 2 \p) 
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The result is then averaged over all configurations of 
the gauge fields. This procedure is repeated for several 
values of the mass m (we have used m — 25 MeV, m = 
50 MeV, m = 75 MeV) , and the obtained dependence on 
m is linearly extrapolated to zero to yield the expectation 
values in the chiral limit. In order to have a well-defined 
limit m — > 0, exact zero modes of the Dirac operator 
were not included in the Dirac propagator. The removal 
of zero modes is legal, since their contribution becomes 
insignificant in the limit V — > oo, m — > if one first 
tends the volume to infinity and then sets the mass to 



zero. The Dirac propagator is evaluated by inverting the 
massive Dirac operator in the subspace spanned by M 
Dirac eigenvectors which correspond to M nonzero Dirac 
eigenvalues with smallest absolute values. The truncated 
operator (J6]) reads as follows: 
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The value of M is limited by the numerical procedure 
(ARPACK in our case) used to find the eigensystem of the 
Dirac operator. We have found that the UV-finite quan- 
tities depend only weakly on M. 

Uniform magnetic field with field strength F\ 2 = B is 
introduced into the Dirac operator by substituting the 
su (2)-valued vector potential A^ with the u (2)-valued 
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The expression @ is valid in infinite space. In order 
to combine it with periodic boundary conditions on the 
lattice we have introduced an additional z-dependent 
boundary twist for fermions, as explained in [25 1. Since 
the total magnetic flux is quantized, on finite periodic 
space with the period L the uniform magnetic field can 
take only discrete values: 
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where q is the quark charge. In our calculations we take 
q = — e/3, corresponding to the d-quark. Due to the 
quantization condition (|10p there is a minimal nonzero 
value of magnetic field, which is equal to y 'qB — 348 MeV 
in most our simulations (for 14 4 lattice with lattice spac- 
ing a = 0.103 fm one has L = 14- a = 1.44 fm). 

In our work the spatial tensor indices are i,j = 1,2,3 
and the Euclidean time is labeled by the index 4. The 
magnetic field F 12 = B 3 = B is uniform in our simula- 
tions, see Eq. ©. 



B. Observables and topology 

In order to characterize the CME quantitatively, we 
studied two basic quantities. The first quantity is the 
local chirality 



P5 (x) = ip (x) 75^ (x) , 



11) 



which is the operator of the difference of the local densi- 
ties of left- and right-handed quarks. 

The second quantity is the electromagnetic current 



(12) 
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All other studied quantities are the powers and the prod- 
ucts of (fTTj) and (H2J). Note that in Euclidean space 
the spinor conjugation is just the complex conjugation 

^Eucl = ^ M- 

According to the orig inal idea of Refs. [l|, |4[ the CME 
appears in the presence of the topologically nontrivial 
configurations. The nonzero topological charge is needed 
to destroy a balance between the left and right quark 
zero modes. In general, configurations of gauge fields 
may have different values of the global topological charge. 
According to the Atyah-Singer theorem, the topological 
charge of a gauge field configuration is given by a differ- 
ence of the right and left zero modes of the Dirac opera- 
tor, 

Nl — N R = Q (13) 

(this is essentially a Euclidean version of Eq. ([2]) written 
for one species of fermions, Nf = 1). The crucial ad- 
vantage of our numerical procedure is that we are using 
the chirally symmetric Dirac operator which possess ex- 
act zero modes. The knowledge of the zero modes for a 
given gluon field allows us to calculate the integer- valued 
topological charge of this configuration using Eq. (TT3")) . 

Surprisingly, we have observed that all studied CME- 
related observables are independent of the global charge 
of the gluon fields. More precisely, in the external mag- 
netic field the effect of the nontrivial topology is smaller 
then the statistical error for any studied quantity like 
(fTTj) and (fT2|) . their powers and correlations. One could 
incorrectly conclude that the electromagnetic properties 
of the QCD vacuum in the external magnetic field have 
nothing to do with the topology of the gluon fields and, 
consequently, the CME is not realized in the (thermal) 
equilibrium (we remind the reader that the CME is an 
interplay between the electromagnetism and the gluon 
topology). 

Coming a bit in advance, we state that the real sit- 
uation turns out to be more complicated. In fact, each 
gauge field configuration contains quite strong local fluc- 
tuations of the topological charge. In the real equilibrium 
configurations of the gauge fields, the density of the topo- 
logical charge is localized on low-dimensional objects [261 ] 
which do not resemble smooth structures like, for exam- 
ple, classical instantons. And, consequently, the effect 
of the global topology is completely overshadowed by 
strong local fluctuations of the topological charge. Thus, 
below we do not classify the gauge field configurations 
by the topological charge, treating all configurations on 
equal footing. On the contrary, our studies, reported in 
this paper, show that the CME does exist in equilibrium 
and that the magnetic fields do affect the electromag- 
netic properties of vacuum at local concentrations of the 
topological charge. 

We prefer to study correlation of the electromagnetic 
current (p~2|) with the chiral density (fTTj) . The chiral den- 
sity is calculated with the help of the low eigenmodes of 
the Dirac operator in the background of the gauge field 
configuration, as described in the previous subsection. In 



fact, the CME predicts the correlation of the electromag- 
netic currents with the imbalanced density of the left and 
right quark zero modes. The exact zero modes of separate 
topological lumps in a configuration of the gluon fields 
(one can think, for example, of separate instantons and 
anti-instantons) become near-zero modes because of the 
interactions between the modes. This interaction breaks 
the zero-mode degeneracy and the former zero modes ac- 
quire tiny corrections. Such near-zero modes carry in- 
formation about unperturbed zero modes including the 
distribution of chirality and, eventually, about the chiral 
magnetic effect. 

In order to illustrate the insignificance of the exact 
Euclidean [27} zero modes let us consider a configuration 
consisting of an instanton-antiinstanton pair. Assume 
that these objects are separated by a distance which is 
much larger compared to their sizes, and this configura- 
tion has no zero fermion modes since the total topolog- 
ical charge is zero. If one applies an external magnetic 
field to this configuration, the lumps of the topological 
charge will be the sources of the generation of the elec- 
tric current. At the instanton the electric current will 
be generated in one direction, while at the antiinstanton 
the current will start flowing in the opposite direction. 
The CME mechanism leads to the current generation be- 
cause each of the (anti)instantons changes the chirality 
of quarks in the positive or negative direction, and the 
local imbalance of the chirality leads to a generation of 
the positive or negative current, respectively. 

At the instanton-antiinstanton configuration the global 
electric current is zero, while locally each topological 
lump induces the electric currents exhibiting the CME. 
Therefore, we concentrate our attention on local corre- 
lators ignoring the global signatures of the CME like 
global currents and/or the Euclidean zero modes. An- 
other argument in favor of local treatment of the CME 
is that in the Euclidean space the exact zero mode pro- 
vides zero contribution to the electric currents and the 
electric charge density. This statement can readily be 
proven in the chiral basis of the Dirac matrices in which 
the 7-matrices are 2 x 2-diagonal structures. Neverthe- 
less, the topological charge fluctuations do induce the 
electric currents because the in the presence on many 
positively and negatively charged lumps of topological 
charge (one can think of, e.g., an ensemble of instantons 
and antiinstantons) the zero modes split into very-near- 
zero modes, which do remember about the local topology 
of the lumps. The low-lying modes do contribute to the 
electric quark's currents and this fact is clearly seen in 
our simulations. 

Thus, in all our calculations we study the nonzero 
fcrmionic eigenmodes and we ignore the exact zero 
modes. Note that in large enough volumes the fermionic 
zero modes are suppressed and in the thermodynamic 
limit their contribution to certain physical observables - 
for example, to the chiral condensate [28j and to the chi- 
ral magnetization [l8| - is zero. Indeed, as the volume 
increases the number of the exact zero modes growing 



as the square root of the total volume while the number 
of the near-zero modes grows proportionally to the first 
power of the volume. 

Finally, we mention that the space-time distribution 
and the dimensionality of the regions of space - where the 
left/right currents are imbalanced - may be as compli- 
cated as the distribution of the topological charge itself. 
The CME is expected to be realized due to the presence 
of the regions where the chiral density (fTTj) of the for- 
mer zero-modes (and now near-zero modes) of quarks is 
enhanced. We can probe the mechanism of the CME 
by investigation of the mutual correlations between the 
chiral density (fTTj) and the electromagnetic current (fT2)> . 

C. Electric charge density: an illustration 

The basic feature of the CME is the spatial charge sep- 
aration, which appears due to the presence of the non- 
trivial magnetic field B. The separation is a result of 
the existence of the quarks' electric current induced in 
the direction of the magnetic field. We call this current 
as the "longitudinal" electric current. According to our 
notations, the longitudinal direction is given by fi = 3 
for T ^ and it covers the whole 03-plane for T = 0. 
The transverse (with respect to the magnetic field) com- 
ponents of the electric current belong to the 12-plane. 

In order to illustrate the CME we have first studied a 
spatial excess of the charge density due to the presence 
of the magnetic field, 

jo (x; A, B) = (jo{x)) A ,B- ( jo (x))a,b=o ■ (14) 

The subtraction in (|14p removes all ultraviolet contribu- 
tions leaving us with infrared nonperturbative contribu- 
tions originating due to the external magnetic field. 

Level surfaces of two typical charge density distribu- 
tions in a fixed timeslice of a zero-temperature gauge field 
configuration are shown in Figures [TJ and for two values 
of the external magnetic field B. 

We can deduce three interesting qualitative features 
from Figures [1] and [5] 

Firstly, we notice that the external magnetic field in- 
deed induces a spatial excess of the electric charge den- 
sity. The stronger the field the larger the excess. This 
property is in agreement with our present understanding 
of the CME. 

Secondly, the excess of the charge density is spatially 
extended along the direction of the magnetic field. This 
feature can be understood in terms of the standard elec- 
trodynamics [2^ |: the quarks are localized in the trans- 
verse plane with respect to the field as they tend to oc- 
cupy low Landau levels (in Figures [TJ and [2] the transverse 
planes are horizontal) . If there were no background gluon 
field then the longitudinal quark trajectories are not af- 
fected by the field. In general quarks have a nonzero 
momentum in the longitudinal direction, so that the tra- 
jectories of the quarks should be elongated strictly in the 
direction of the magnetic field. The gluon background 
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q B = 0.7 GeV 2 




FIG. 1: (color online). The excess of the density of the electric 
charge (|14p due to the applied magnetic field qB = 0.7 GeV 2 
on 14 4 lattice with the volume (1.44 fm) 4 . A typical three di- 
mensional slice of a typical gauge field configuration is shown. 
The magnetic field is directed vertically from the bottom of 
the picture to the top. The red regions mark the excess of 
the positive charges while the blue regions correspond to the 
negative electric charge density. 



qB = 1.8 GeV 2 




FIG. 2: (color online). The same as in Figure (JTJ) but for qB = 
1.8 GeV 2 and for the configuration of non-Abelian gauge field. 



field provides a visible distortion to the quarks' motion. 
This feature is observed in both figures. 

Thirdly, Figures [TJ and [2] illustrate a conceptual diffi- 
culty associated with identifying the electric dipole mo- 
ment of the spatially separated charges in a finite volume. 
Due to the (periodic) boundary conditions the electric 
dipole moment is not a universal notion as it depends 
on the choice of the reference frame. Consider, for ex- 
ample, one dimensional space of the finite length L with 
periodic boundary condition imposed. Two points sep- 
arated by the distance r > have the dipole moment 
fj, e = er > 0. However, if we change the direction of our 
reference frame, one can count the distance between the 
two points as L—r. As a consequence, we get the different 
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value for the electric dipole moment p, e — e(r — L) < 0. 
The difference between the two definitions results in the 
large uncertainty of the electric dipole moment, and this 
fact makes the physical sense of the definition obscure. 
One can, of course, limit all distances on the lattice by 
r ^ L/2, so that if r < L/2, then we define pb e = er 
while r > L/2 then p e = e(r — L). However, this defini- 
tion is not satisfactory as well, because small variations 
of the charge around the most distant point r ~ L/2 
lead to large variations of the electric dipole moment, 
fi e ~ eL/2 — > —eL/2. Another solution of this problem 
could be to replace the infinite- volume definition p e = er 
by its finite- volume counterpart, ji e = ef(r, L), where the 
function /: (i) obeys the symmetry f(r, L) — f(L — r, L) 
and (ii) reduces to / — > r in the infinite volume limit, 
L — > oo. However, this solution is not satisfactory ei- 
ther, because the electric dipole moments in intermedi- 
ate volumes would inevitably be dependent on the form 
of chosen function f(r, L). 

Thus, the electric dipole moment is an ambiguous no- 
tion in finite volume. Therefore, instead of the electric 
dipole moments we study the electric currents. The in- 
duced currents represent a complimentary aspect of CME 
since the emergence of the electric dipole moment is nec- 
essarily associated with an electric current The advantage 
of the current-based approach is clear: the definition of 
the electric current (TT2|) is a local notion which is not 
subjected - at least directly - to the finite volume ambi- 
guities. 



D. Chirality fluctuations and CME in cold matter 

Since the CME originates due to fluctuations of the 
local chirality, we first study the effect of the external 
magnetic field on the expectation value of the square of 
the local chirality (JTTJ) : 



(pt) 7B (B,T) = i J d 4 x( P5 (x) P5 (a) 



B,T 



-yjd i x{p b (x)p b {x)) B =o,T=a- (15) 
v 

Here V is the total volume of the lattice, and (...) b.t 
means expectation value with respect to (w.r.t.) the 
thermal state at temperature T with background mag- 
netic field B. Subtraction of the expectation value at 
zero temperature (T = 0) and in the absence of the ex- 
ternal field (B — 0) removes all ultraviolet divergences 
providing us with a physical quantity which should be 
insensitive to the ultraviolet cutoff. In fact, the quan- 
tity ([15)) provides us with the strength of the chirality 
fluctuations because in the equilibrium QCD the average 
of the chirality is zero, (p 5 ) = 0. We denote the result of 
the subtraction by the subscript "IR" in ( . . . ) //? stress- 
ing the importance of the nonperturbative infrared (IR) 
contributions. In all our expectation values the corre- 
sponding zero-temperature and zero-field contributions 



are always subtracted as it is done in (|15p . Note that 
the coordinates of all operators entering our observables, 
including Eq. (fT5"|) , are taken at one point x. 
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FIG. 3: (color online). The expectation value of the chirality 
squared (|15[) . pi, vs. magnetic field at different lattice param- 
eters (the lattice spacing a and the lattice volume L 4 ), and 
for different number M of Dirac eigenmodes which were taken 
to determine the truncated Dirac propagator ([8]) used in Q 
to calculate (|15p . 

Due to the subtraction, the expression (fTS")) gives us 
the influence of the magnetic field and/or finite temper- 
ature on the strength of the chirality fluctuations. We 
plot this quantity in Figure [3] for a few values of the 
lattice spacing a, lattice volumes L 4 and the number of 
low eigenmodes M which are used to calculate the Dirac 
propagator ([5]) . We found that the variations of all these 
parameters change the expectation values by less than 
10%, thus confirming good scaling of our results. 

The CME appears due to the local topology of the 
gauge field configurations. The topological content of 
the gauge field configurations is reflected in the struc- 
ture of the exactly-zero and near-zero eigenmodes. The 
exactly zero eigenmodes are not significant in the Eu- 
clidean space as we have discussed above. However, the 
near-zero modes play an important role in the CME, and 
they also come with larger weights in the propagator l[5]). 
Thus, the consideration of only low-lying zero modes is 
justified. For our particular observables, the errors asso- 
ciated with our truncation scheme are smaller than the 
typical statistical errors. However, we should make a 
cautionary remark that as the volume of the lattice in- 
creases, the density of low-lying eigenmodes grows, and 
one should take into account more and more of them. 

Notice that most of Figures in our paper are plotted in 
the logarithmic scale, and that due to dimensional rea- 
sons the data for chirality fluctuations and, later, for cur- 
rent fluctuations are always expressed in GeV 6 . Due to 
the high (sixth) power the actual values of our quantities 
look quite small, but this fact in no way indicates that 
the physical effects are small. For example, in Figure [3] 
at qB = 1.5 GeV 2 the value of the chirality fluctuation is 
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FIG. 4: (color online). The expectation values of the chi- 
rality squared (|15|l vs. the magnetic field at three different 
temperatures. 



(/?§) = 10~ 4 GeV 6 . This number corresponds to the large 
quantity w (215 Me V) 6 if the magnitude is expressed 
in terms of MeV. 

In Figure 2] we show the dependence of the expecta- 
tion value (p\)xr{B,T) on the magnetic field at three 
temperatures: at zero temperature, at finite tempera- 
tures in the confinement phase (T = 0.82 T c ) and in 
the deconfinement phase just above the phase transition, 
T = 1.12 T c . At zero temperature the fluctuations of the 
chirality quickly grow with magnetic field. As the tem- 
perature increases the grows rate gets smaller and in the 
deconfinement phase the chirality fluctuation are almost 
insensitive to the value of the magnetic field. 

In the absence of the external magnetic field, B = 0, 
the rise of temperature leads, according to Figure |H to 
the increase of the chirality fluctuations. However, suf- 
ficiently strong magnetic fields can also induce the en- 
hancement effect comparable to the effect of the thermal 
fluctuations. For example, the chirality fluctuations at 
zero temperature in the magnetic field of the strength 
qB ~ 2 GeV 2 has the same value as the fluctuations of 
chirality in the absence of the external magnetic field at 
temperature T = 1.12 T c : 



(/0§(x)) T=0) B= 2GcV 2 ~ (pI{x))t=1. 



12 T c , B=0 



(16) 



It is known that the temperature suppresses the fluc- 
tuations of the global topological charge [30], 



xq = ^(Q 2 ) 



(17) 



Lattice studies show that in Yang-Mills theory the fluc- 
tuations are almost independent on the temperature in 
the confinement phase, 

Xq (T = T c )^ X q(T = 0), 

while as temperature exceeds the critical point, T > T c , 
the fluctuations vanish quickly [3lf | . For example, in 



SU(2) Yang-Mills theory just above the phase transition, 
T = 1.1 T c , the value of the fluctuations is more than 
twice smaller compared to the value in the confinement 
phase [311 ] . 



x Q (T = i.ir c )«o.5x Q (r = T c ). 

Moreover, it is clear that in quenched theory the topo- 
logical charge of the gluon fields is not affected by the ex- 
ternal magnetic field, because the coupling between glu- 
ons and electromagnetism is provided by the dynamical 
quarks which are absent in the quenched theory. These 
observations suggest that the local fluctuations of the chi- 
rality (|T5|) and the fluctuations of the global topological 
charge (jTTJ) are not related directly. It is more likely 
that at high temperatures chirality-changing transitions 
occur due to local processes involving sphaleron-like con- 
figurations 32]. At zero temperature one might explain 
chirality changes as follows: in a strong magnetic field 
quark undergoes one-dimensional motion along the mag- 
netic field. Quark chirality can be efficiently changed in 
local one-dimensional scattering processes on the fluctu- 
ations of gauge fields, where the spin is not affected but 
the momentum changes its sign. In the deconfinement 
phase this mechanism does not work because the influ- 
ence of the gluon fluctuations on the quark's motion is 
lower in this phase (the quark confinement is lost, for 
example) and one can suggest that the typical scattering 
does not reflect the momenta. 

Concluding this subsection, we notice that the CME 
crucially depends on the imbalance between the left and 
right chiral modes. The local strength of the chiral im- 
balance is characterized by the value of the chirality fluc- 
tuations, (p 2 ,), which determines the rate of chirality- 
changing transitions. Since the external magnetic field 
strongly enhances the chirality fluctuations (they become 
comparable with the fluctuations in the deconfinement 
phase) so that the Chiral Magnetic Effect - if we imagine 
for a moment that the property of the color confinement 
were either absent or inessential due to density effects - 
could be observed in cold nuclear matter as well. 



E. Fluctuations of electromagnetic current 

We checked that the expectation value of all compo- 
nents of the electromagnetic current itself (fT4")) is zero 
within the error bars. This fact is not unexpected be- 
cause the magnetic field passes through the regions of 
space-time characterized by positive and negative chirali- 
ties (fm . According to the CME, in regions with opposite 
chiralities the magnetic field induces the electric currents 
in opposite directions with respect to the magnetic field. 
Thus, globally, the longitudinal electric currents should 
cancel each other. As for the transverse (with respect to 
the magnetic field) spatial directions, the electric currents 
should cancel each other because of circular motion on 
(perturbed) Landau levels. Finally, the average timclike 
current is zero because of the Gauss theorem (the total 
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FIG. 5: (color online). The same as in Figure [3] but for the 
fluctuations (fT51) of the electromagnetic current j'3 (fT2|) . 

electric charge in a closed volume with periodic boundary 
conditions should be zero). These statements are valid 
both at zero and at non-zero temperatures. 

Below we study the dispersion of the electromagnetic 
current 

Ul)m(B,T) = ± J d^jl (x))b,t 
v 

-y J d A x{jl (x) ) b =o,t=o [no sum over fjt] , (18) 
V 

instead of the expectation value of the electric current. 
Note that here and below no sum over the Lorenz index 
/1 is assumed in all our expressions involving . 

Following study of the fluctuations of chirality, we first 
check that the internal lattice parameters do not influ- 
ence the final results in Figured Then in Figure[S]we plot 
the expectation values of the fluctuations of each compo- 
nent of the electric current at zero temperature. At T = 
the rotational symmetry is broken only by the electro- 
magnetic field strength tensor with F\% 7^ 0, so that 
the components of the electric currents can be grouped 
into transverse and longitudinal sets, respectively, 

Ul) = {fi), (io 2 ) = (ii) (T = o). (19) 

At weak magnetic fields all components have the fluc- 
tuations of the same order. However, at stronger fields 
the fluctuations of the longitudinal currents (j'3 and jo) 
are stronger compared to the fluctuations the transverse 
currents (jx and j-i). Both transverse and longitudinal 
components grow with magnetic field, but the longitu- 
dinal components grow faster. Thus we see that at zero 
temperature the magnetic field enhances both the current 
j 3 in the direction of the magnetic field and the charge 
density jo, which is a clear signature of the charge sep- 
aration associated with the CME. The enhancement of 
the transverse components j\ and j% in magnetic field 
is natural since the transverse momentum of a particle 



FIG. 6: (color online). The expectation value of the squares 
of the transverse (ji, j'2) and longitudinal (j'3, jo) components 
of the vector current ()12[) vs. the applied magnetic field at 
zero temperature. 
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FIG. 7: (color online). The same as in Figure [6] but for T = 
0.82 T c . The longitudinal components j$ and jo differ due to 
nonzero temperature (|20[1 . 

occupying a Landau level grows with the increase of the 
magnetic field. 

At nonzero temperature the (discrete) rotational sym- 
metry is broken not only by the electromagnetic field, 
F12 =/= 0, but also by the compactified temperature di- 
rection fi = 0. Therefore the charge fluctuations {ji — 
0) loose the similarity with longitudinal (/i = 3) cur- 
rent fluctuations, while the fluctuations of the transverse 
(fx = 1, 2) currents are still degenerate, 

(ii 2 ) = (ii), (j 2 )^(ji) (7V0). (20) 

In Figure [7] we plot the fluctuations of all components 
of the electric currents at finite temperature T — 0.82 T c , 
at which the system is still in the confinement phase. One 
can see that at zero magnetic field the currents fluctuate 
strongly due to thermal fluctuations. All components 
of the electric currents fluctuate approximately with the 
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FIG. 9: (color online). The ratio of the longitudinal current 
to the transverse one vs. the strength of the magnetic field at 
various temperatures. 



same strength. As the magnetic field grows the fluctua- 
tions of all components drop a bit, and then they start 
to grow steadily. The grows rate is noticeably smaller 
compared to the slope observed at zero temperature, Fig- 
ure [SJ On the other hand, there are similarities with the 
zero-temperature case: (i) the fluctuations of the charge 
density and of the longitudinal (with respect to the mag- 
netic field) current have equal magnitudes; and (ii) the 
longitudinal currents fluctuate stronger compared to the 
transverse currents. 

The situation turns to be quite different in the decon- 
fmement phase, Figure [H First, the fluctuations of the 
longitudinal current are insensitive to the strength of the 
magnetic field. Second, both the charge and transverse 
current fluctuations are decreasing functions of the field 
strength. The drop in the electric charge fluctuations is 
especially noticeable. This feature can be explained by 
an effect of the thermal Debye screening [33| . 

In Figured we show the ratio of the longitudinal com- 
ponent of the current to its transverse component. At 
zero temperature the current is substantially enhanced 
at qB ~ 1.5 GeV 2 , while at relatively weak fields, qB ~ 
0.3 GeV 2 , the enhancement of the current is not seen. 

At nonzero temperature the relative enhancement is 
weaker. The ratios at T = 0.82 T c and at T = 1.12 T c 
have similar properties which are in a qualitative agree- 
ment with the CME features: at low magnetic fields the 
longitudinal and transverse directions are equivalent and 
the increase of the magnetic field gives a bit stronger 
preference for the longitudinal direction with respect to 
the transverse one. 

We have already mentioned the reasons why the CME 
is expected to exist in the high temperature phase [l|, 
where the confinement is absent and the chiral symmetry 
is restored. The reasons are quite simple: the presence 
of the chiral condensate destroys imbalance between the 
right and left chiral modes while the confinement does 
not allow the spatial separation of the quarks at the dis- 
tances larger than the confinement distance. This ar- 



gument seems to be independent on the inclusion of the 
dynamical (virtual) quarks into the consideration, so that 
it should be valid, in lattice language, both in quenched 
and in unquenched theories. 

We have found that in the quenched approximation the 
qualitative signatures of the chiral magnetic effect exist 
both in the confinement and in the deconfincment phases. 
In all our simulations, at zero temperature and at finite 
temperature above and below the critical point, we see 
that the longitudinal component of the electric (vector) 
current dominates over the transverse component. How- 
ever, in the confinement phase the electric current rises as 
the function of the external magnetic field while in the 
deconfinement phase the value of the current is almost 
insensitive to the strength of the magnetic field. 

F. Chirality-current correlations 

According to the CME the longitudinal electric cur- 
rent may locally be induced by the imbalance of chirality. 
Therefore it is natural to expect that the local fluctua- 
tions of the longitudinal electric currents are correlated 
with the chirality fluctuations. However, the correlator 
of the chirality and the electric current ( ps ) vanishes, 
since in isotropic homogeneous space it is impossible to 
construct an axial vector from an antisymmetric tensor 
of magnetic field strength. Therefore, we have consid- 
ered the correlation between the squares of the induced 
current and the chirality: 

c(P 5 ,J»)~ {4W) • (21) 

It turns out that in the confinement phase the corre- 
sponding correlation functions for all components of the 
electric currents vanish within statistical uncertainty, 

c(pIjI)c0 T<T c . (22) 
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In the deconfincmcnt phase the correlation function (|2"T|) 
is not zero, Figure QJJ while the correlation decreases 
with the increase of the magnetic field. All these facts 
suggest that the enhancement of the fluctuations of cur- 
rent in the direction of the magnetic field is not locally 
correlated with the chirality fluctuations. 

We would like to emphasize, that the absence of the lo- 
cal correlation does not in general mean the independence 
of the longitudinal currents on the chiral density. The 
longitudinal current can be caused by the chirality fluc- 
tuations, and, at the same time, the current itself may be 
not localized in the local regions of the spacetime where 
the chiral density is enhanced. Although this statement 
sounds puzzling, similar effects are already known in the 
condensed matter physics. For example, the well-known 
Anderson localization [34[ of electron inside a semicon- 
ductor is caused by the presence of impurities. However, 
the spatial regions, in which the electron wavefunctions 
are localized, do not coincide with the positions of the 
impurities themselves. 



G. Clean case: instanton-like configuration 

In the previous sections we have studied the CME us- 
ing the gluon fields of the quantum vacuum. Such gluonic 
fields contain ultraviolet fluctuations both of perturba- 
tive and nonperturbative nature. In particular, the topo- 
logical charge density is subjected to strong fluctuations 
with quite nontrivial localization properties, so that the 
CME should be realized locally in a vicinity of each topo- 
logical lump. In order to illustrate the existence of the 
CME in a more theoretical fashion, we have prepared a 
smooth instanton-like gluon field with a unit topological 
charge on the lattice 14 4 and applied the external mag- 
netic field to this configuration. In Figure [TJJ we show 
the behavior of the chirality fluctuations and the fluctua- 
tions of longitudinal electric current at this instanton-like 



FIG. 11: (color online). The fluctuations of the chirality (|15[) 
and the fluctuations of the longitudinal electric current (|18p 
at an instanton-like configuration with unit topological charge 
as a function of the applied magnetic field. 

background. 

The first indication of the existence of the CME comes 
from the observation of the increase of chirality fluctua- 
tions with increase the strength of the external magnetic 
field. Figure [TJJ illustrates that the external magnetic 
field, on average, enhances a local imbalance between 
approximately left and approximately right fermionic 
modes in the background of this Q — 1 configuration. 
The effect may readily be understood as a consequence 
of the CME: the magnetic field tend to localize "left" 
and "right" fermionic modes in different parts of the 
space-time because the left and right modes have op- 
posite magnetic moments. The external magnetic field 
acts as a chemical potential which provides a different 
weight to these modes. The left and right modes have 
opposite direction of the momentum with respect to the 
magnetic moment so that the external magnetic field tries 
to separate them spatially as well. Thus, the positive and 
negative chiralities do not compensate themselves in p§, 
leading to increase of the chirality fluctuations. 

The second signature in favor of the CME appears in 
the properties of the (squared) longitudinal electric cur- 
rent. Notice that the "longitudinal" current contains, in 
fact, contribution from the two components, j'3 and jo, 
because they are indistinguishable from a geometrical 
point of view (we remind that we are working at zero 
temperature in external magnetic field B = F^). There- 
fore we also plot in Figure [TT] the absolute value of the 
longitudinal current squared, jjj = jg + jf which turns 
out to be an increasing function of the magnetic field. 
This behavior is consistent with our understanding of 
the CME. 

In order to provide an illustration of the effect, we vi- 
sualize the squared transverse and longitudinal electric 
current densities for the Q — 1 configuration in the ex- 
ternal magnetic field. In Figures [T2] and [T3] we show the 
current densities in the 12 and 34 cross-sections, respec- 
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FIG. 12: (color online). The squared components of the elec- 
tric current (in arbitrary units) in a 12-plane on the back- 
ground of an instanton-like configuration with topological 
charge Q — 1. The upper (green) sheet represents the spa- 
tial distribution of the longitudinal current jf, = jo + j'f , and 
the lower (violet) sheet corresponds to the transverse current 
3x =3i +J2- 




FIG. 13: (color online). The same as in Figure [121 but for a 
03-plane of the same gluon field. 



tively. One can clearly see that in a strong magnetic 
field the longitudinal component is always larger than 
the transverse one in agreement with the prediction of 
the CME. 

Concluding this subsection we summarize, that our 
study of an instanton-like configuration of the gluonic 
field with the topological charge Q = 1 supports the ex- 
istence of the Chiral Magnetic Effect. 



particles produced in each event, and <pi a , <f>jb are the az- 
imuthal angles of the momenta of the produced particles 
w.r.t. the reaction plane. 

These observables can be expressed in terms of the 
total charge separated in each event as: [l|: 



c(AaA b ) 
(N a )(N b ) 



(24) 



where a, b — ±, A + is the difference of the total numbers 
of positively charged particles emitted above or below 
the reaction plane, A_ is the same quantity for neg- 
atively charged particles and c is some coefficient (of 
order of unity) which depends on the characteristics of 
the hydrodynamical flow in the reaction. One also has 
( N a ) = (Nb) = N q , where the average total number of 
the particles of the same charge per event N„ (for Au-Au 
collisions) can be found, e.g., in Table 1 in [lj. 

In order to compare experimental data with our re- 
sults, it is convenient to consider the following combina- 
tion of a a b'- 



a ++ + a — 2a_| 



(AQ) 



(A + -A_) 



,(25) 



where AQ = A + — A_ - is the difference of total charges 
of particles emitted above and below the reaction plane. 
( (AQ) ) can be estimated from our lattice data on ( jjj ) 
using the following rough model: consider a fireball of 
radius R created during the heavy-ion collision, and as- 
sume that it emits positively or negatively charged par- 
ticles from each element of surface <5<7j, with \8<Ji\ ~ p 2 . 
The quantity p is some characteristic correlation length, 
so that emission of particles on each surface element of 
the area p 2 is statistically independent. AQ is calculated 
as the difference of the total charges emitted in some 
characteristic collision time r from the upper and from 
the lower hemispheres of the surface of the fireball: 



SQ 



8a \ ■ B I ji- 8ff l 



(26) 



For the expectation value ( (AQ) 2 ) one then has: 



((AQ) 2 )=T 2 J2(U-5* 



(27) 



III. CHARGE ASYMMETRY: COMPARISON 
WITH EXPERIMENT 



rV£<ft^) )\Sa t \=r 2 p 2 / dn{j-n) 



Experimental signatures of the Chiral Magnetic Effect 
are usually discussed in terms of the following observ- 
ables, suggested first in Ref. 0: 



^ = fEf«EE cosI 



, N aNb ■ , , 



+ <)>jb) , (23) 



where a, b = ± labels positively/negatively charged par- 
ticles, N e is the total number of experimental events, N a 
and Nb are the numbers of positively/negatively charged 



We further assume that emission in different directions 
is independent: (j^jv) = for /i / c. With these as- 
sumptions, we finally arrive at the following expression 
for((AQ) 2 ): 



((AQ) 2 
N 2 



4lTT 2 p 2 R 2 

3N 2 



J\\) 



'A 



(28) 



where j'ii and j± are, respectively, the currents along the 
magnetic field and perpendicular to it. We take the val- 
ues of N q from Table 1 in [1[ and use our lattice data on 
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FIG. 14: (color online). The average fluctuations of the dif- 
ference of total electric charges of particles emitted above and 
below the reaction plane: comparison of the experimental ob- 
servable of 0, 0| (the filled triangles) with our lattice data 
(the filled circles). 



(j 2 ) and ) at T = \.12T C . Here we take into account 
the fact that the expectation value of the squared spatial 
currents in the Minkowski and Euclidean spaces are the 
same. 

In Table 1 in |l| one can also find a relation between 
centrality ("% most central") and the ratio b/R, where 
b is the impact parameter. From expression (A. 12) in 
[l| one can roughly estimate the magnetic field during 
the collision as eB ~ (0.16/i?) GeV 2 . For simplicity, we 
assume that the magnetic field is nonzero and constant 
only during the time r. In order to obtain ( (AQ) 2 ) at 
small values of magnetic field we have used polynomial 
interpolation. The resulting plot is shown in Figure 1141 
Experimental points obtained by STAR collaboration are 
taken from the plots presented in @, S] • The coefficient 
t 2 p 2 R 2 in (|28| was fitted so that the lattice data and the 
theoretical data match best. One can take, for example, 
t ~ 1 fm, R ~ 5fm, p ~ 0.2 fm, which are reasonable 
values by order of magnitude, especially for such a rough 
model. 

We observe a good agreement between our data and 
the experimental results. Despite the agreement is seen 
at the quantitative level, the approximations used in our 
lattice technique suggest that this result should be consid- 
ered as a qualitative indication rather than a solid proof. 
It should be noted that at T — 1.1 2T C the dependence 
of ( j 2 ) and hence of ( (AQ) 2 ) on the magnetic field is 
rather weak, and, moreover, the current does not increase 
with magnetic field. We have to conclude that the dom- 
inant effect which determines the shape of the plot on 
(fTH) is the decrease of N q with centrality (see Table 1 
in [l[). It could be therefore much more informative to 
consider the observables b a b introduced in [l|, which do 



not contain the normalization factors 



N a N b ■ 



These ob- 



servables can be directly compared with our data for the 
observable ( j 2 ) , both in the deconfinement phase and 



at zero temperature. Note also that since at zero tem- 
perature the current increases with magnetic field much 
faster than at high temperature, it might be not unrea- 
sonable to search for strong magnetic effects also in cold 
(and dense) nuclear matter. 



IV. CONCLUSIONS 



In this paper we studied possible signatures of the Chi- 
ral Magnetic Effect [E IIS, El in quenched SU (2) lattice 
gauge theory using the chirally invariant Dirac opera- 
tor. We observed the expected enhancement of the elec- 
tric current in the direction of the magnetic field both 
in quantum configurations of the gluon fields and in a 
fixed instanton-like gauge fields with nonzero topological 
charge. The CME at the instanton-like configurations is 
much more pronounced compared to the real gluonic vac- 
uum (an illustration can be found in Figures [T2l and [T5|) . 
The CME is based on the imbalance between left and 
right fermionic modes. The measure of the local imbal- 
ance is the chiral density of the quarks. On average the 
global chiral charge is zero, while the local fluctuations 
of the chiral density can be quite strong. We found (Fig- 
ure 2]) that at zero temperature the chiral fluctuations 
are drastically (by two orders of magnitude) enhanced in 
the applied magnetic field. 

As the temperature increases the chirality fluctuations 
increase and, simultaneously, they become less sensitive 
to the strength of the magnetic field. In the deconfine- 
ment phase the local chirality fluctuates much stronger 
compared to the fluctuations at zero temperature while 
the fluctuations themselves are practically independent 
on the strength of the magnetic field. The enhancement 
of the chirality fluctuations due to magnetic field can be 
illustrated by the following example: in strong magnetic 
field qB = 2.5 GeV 2 at zero temperature T = the chi- 
rality fluctuations take the same value (p 2 ) « (220 MeV) 6 
as in the case of zero magnetic field, qB = and high 
temperature, T = \.Y1T C . 

In the confinement and in the deconfinement phases 
the properties of the induced electric currents differ 
strongly. At zero temperature the longitudinal (i.e., di- 
rected along the magnetic field) currents grow with the 
increase of the strength of the magnetic field. This ef- 
fect is accompanied by a (weaker) enhancement of the 
transverse currents, which we attribute to a transverse 
squeezing of the Landau levels (Figure . 

At nonzero temperature but still below the transition, 
T = 0.82T C , the slopes of enhancement of both longi- 
tudinal and transverse electric currents become smaller. 
However, the longitudinal current is enhanced stronger 
than the transverse current (Figure [7|). 

In the deconfinement phase at T = 1.1 2T C the fluctua- 
tions of the longitudinal currents are almost independent 
on the magnetic field. Moreover, the transverse compo- 
nents of the electric current decrease slightly as the mag- 
netic field gets stronger. The fluctuations of the electric 
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charge density become drastically suppressed in the de- 
confinement phase in the strong enough magnetic field 
(Figure IE). 

We also found that the presence of a nonzero global 
topological charge in the real gluonic configurations leads 
to a little effect on the longitudinal currents. This hap- 
pens because the influence of a global charge is largely 
overshadowed by strong local fluctuations of the topolog- 
ical charge. We conclude, that in our simulations the 
CME is realized at local lumps of topological charge in 
the gluonic configuration. 

We would like to make here a cautionary notice that 
our investigation utilized many approximations. We used 
two colors instead of three, our simulations are done in 
the quenched approximation, the fermionic propagators 
were truncated, and the system was studied in a finite 
volume. We found that the systematic errors associated 
with truncation, finite-size and finite-volume effects are 
smaller than our statistical errors, and we expect on gen- 
eral grounds that the reduced number of colors should not 
affect our results drastically. However the absence of the 
dynamical fermions may in general correct our results. 
We got a qualitative indication of the existence of the 
CME, while the quantitative features of the CME may 
be accessible only in simulations with dynamical quarks. 

Finally, in real heavy ion collisions the magnetic field 
is perpendicular to the reaction plane. Consequently, the 
CME implies an imbalance between total electric charge 



observed above and below the reaction plane which is 
indeed seen in preliminary data published by STAR col- 
laboration at RHIC 0, HI- Using a simple model of a 
fireball we obtain a good agreement between our data 
and experimental results of STAR Collaboration (Fig- 
ure HH) on the average fluctuations of the difference of 
total charges of particles emitted above and below the 
reaction plane. 
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